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As a new scheme of treating the tensor interaction of the nucleon-nucleon interaction, there is a 
proposal of a tensor-optimized shell-model (TOSM) for the study of medium and heavy nuclei. The 
TOSM includes the deuteron-like tensor correlation and provides quite a good reproduction of the 
binding energy and the size of ''He as compared with rigorous few body calculations. We propose a 
tensor-optimized few-body model (TOFM) using the spirit of the TOSM approximation in the few 
body framework with bare nucleon-nucleon interaction. We find that the TOFM can account for 
the strength of the tensor interaction very well and almost reproduces the full energy and various 
energy components as compared with rigorous few body calculations. 

PACS numbers: 21.45.-v,21.10.Dr 



It is very important to describe nuclear structure by 
using bare nucleon-nucleon (NN) interaction. Such an 
effort has been made for a few body [Ij] and light mass 
nuclei up to about A~12 0. The difficulty of solv- 
ing nuclear ground states by using the NN interaction 
is the presence of the strong short range repulsion and 
the medium range tensor interaction caused by the pion 
exchange. The calculated results on the deuteron indi- 
cate that the s-wave component has a distinct dip at the 
short distance and the probability of the d-wave compo- 
nent is of order of 5%. Although the d-wave component is 
small, the dominant attraction of order of 80% is caused 
by the tensor interaction to provide the small binding 
energy by canceling with a large kinetic energy for the 
case of the AV8' interaction. These features have to be 
handled in any theoretical frameworks for a quantitative 
account of nuclear states. The calculated results of few 
body frameworks provide extremely a good account of 
nuclear ground states and a few excited states by includ- 
ing a phenomenological three body interaction 14I . 

In recent years, two very important methods were pro- 
posed for the description of finite nuclei. One is the uni- 
tary correlation operator method (UCOM) and the other 
is the tensor optimized shell model (TOSM). The UCOM 
is a method to handle the short range correlation in terms 
of the unitary correlation operator by truncatingthe re- 
sulting many body operators at the 2 body level [3] . The 
TOSM is a method to treat the strong tensor interaction 
in the 2p-2h model space in the shell model basis [1]. 
Hence, we are able to calculate medium and heavy nu- 
clei by using bare NN interaction. We discuss here the 
essence of the TOSM approximation. The strong tensor 
interaction acting on two nucleons in a spin-saturated 
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shell model state oughts to excite them into two parti- 
cle states. Hence, it is a minimum requirement to take 
2p-2h excitations in order to treat the tensor interaction, 
which provides a large attraction to form nucleus. An- 
other important feature is to use various range (short 
range in particular) gaussian wave functions so that the 
variational function is able to describe the strong tensor 
correlation. Myo et al. introduced the TOSM in order 
to describe the spin-orbit splitting effect in ^He due to 
the blocking mechanism of the Pxl2 orbit caused by the 
deuteron-like tensor correlation [5|. Another important 
application of the TOSM was the formation of the halo 
structure in ^^Li Q. Again the tensor correlation is re- 
sponsible to disfavor a configuration of 2 neutrons in the 
Pi/2 orbit and makes the contribution of 2 neutrons in the 
Si/2 orbit large to be the source of the halo structure. 

A theoretical study was then performed by Myo et al. 
to use a bare NN interaction in the TOSM for ''He and the 
numerical results were compared with those of few body 
methods They used the central correlation part of 
the UCOM in order to treat the short range repulsion in 
the NN interaction. The comparison was quantitatively 
very successful. The tensor matrix element comes out to 
be about 80% of that of few body methods. The binding 
energy was about 3 MeV smaller than that of the few 
body methods. Hence, the TOSM can be used for the 
discussion of nuclear structure of medium and heavy nu- 
clei, in particular for the study of the effect of the tensor 
interaction for various observables in those nuclei. It is 
then very important to know from where the difference 
between the TOSM and the rigorous calculations comes 
from and what should be done in order to improve the 
TOSM description of nuclei. 

In a few body systems, we have a powerful few-body 
technique to describe light nuclei by using the rela- 
tive coordinates, which correspond to the coordinate of 
the nucleon-nucleon interaction. It is much easier to 
take into account the tensor interaction and the short 
range repulsion. On the other hand, the necessary anti- 
symmetrization of all the particle coordinates prohibits 
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the appUcation of the few-body method to larger mass 
nuclei. One of the reasons of the difficulty is a large vari- 
ational model space in which the whole energy is min- 
imized. Hence, it may be a good idea to introduce the 
spirit of the TOSM approximation in the few body frame- 
work. This idea corresponds to take the s-wave configu- 
ration as the basis and add a single d-wave in a relative 
coordinate and do not take double and/or triple d-wave 
configurations. The purpose of the present study is then 
to develop a tensor-optimized few-body model using the 
spirit of the TOSM and compare with results of rigorous 
calculations in a few body framework and to obtain a 
good approximate method to study larger mass nuclei. 

We state first the tensor optimized shell model 
(TOSM), which was developed for the description of the 
tensor interaction in the shell model basis. The TOSM 
wave function is written as 

\'f)^Co\Q)+J2^c.\2p-2h:a) . (1) 

a 

Here, |0) denotes the Op-Oh shell model state. The 2p- 
2h states \2p — 2h : a) are excited by the operation of 
the tensor operator Si2{x) oc [12(2:) x [cti x tT2]2]o on two 
particles in the shell model Op-Oh state. In the tensor op- 
erator, X in the spherical harmonics function Y2 denotes 
an angle of relative coordinate, and cr's are the spin op- 
erators of two particles. Hence, the two components are 
connected by the tensor operator as 

(2p-2/i:a|5i2(x)|0) ^0 , (2) 

where the quantum number a denotes all the possible 
2p-2h states connected by the tensor operator. The 2p- 
2h states are constructed to treat the tensor correlation 
by taking various gaussian range wave functions, partic- 
ularly those with small size to optimize the role of the 
tensor interaction. Hence, the essential approximation 
of the TOSM is the inclusion of the variational states, 
which are excited by only one operation of the tensor in- 
teraction. As for the short range correlation, they use the 
central correlation part of the UCOM ^ in the TOSM 
framework, because it is difficult to express the short 
range behavior in the shell model basis [J]. 

On the other hand, a few body method uses the rela- 
tive coordinates and hence it is more efficient to express 
the tensor and short range correlations, since the tensor 
interaction and the short range repulsion are expressed 
by the relative coordinates. We only have to introduce 
all the states connected by one operation of the tensor in- 
teraction and take enough gaussian functions to express 
both the tensor and short range correlations. We should 
end up introducing the following tensor-optimized wave 
function for A = 2, 3 and 4 body systems by taking the 
essential features of rigorous calculations in the few body 
method. Hence, we write a few body wave function as a 
linear combination of S and D wave components, 

I*) = + \^d) . (3) 



The condition of the D-wave component should be 

(*d|^i2(x)|*s) 7^0 . (4) 

Hence, |^_d) should contain l2-function only once. Par- 
ticularly the most essential D-wave state is to introduce 
a y2-function in the Jacobi coordinate Xi = ri — r2 and 
perform all the necessary permutations so that all the 
particle pairs profit the use of the strong tensor interac- 
tion. 

The S-wave component for s-shell nuclei with A < 4 is 
written as 

Ns 

= E CiAlM'Bs, {xr})Xsi{a^})xti{r.})) ■ (5) 

The spin wave function with total spin s for a A— 4 sys- 
tem is written as 

XsiW;}) = [[[Xi(f^l) X Xi(cr2)]si2 X Xi('73)]si23 (6) 
XXi(c^4)]s . 

The iso-spin wave function is written similarly as 

Xtiin}) = [[[Xi(n) X xi(t2)]u, X xi(r3)]ti,3 (7) 

XXi(T4)]t . 

On the other hand, the spatial wave function is written 
as (1 

M'Bs, {x^}) = exp(-ixBsx) (8) 

Here, x represents the relative coordinate vector x = 
(xi,X2, ..) and xBgx means the short-hand notation of 
^ B^s^i • Xj. The relative coordinate vectors x^ 
with i = 1, .., Np — l are expressed in terms of the particle 
coordinates with j = 1, .., Np, where Np is the number 
of particles of a nucleus. This variational function is able 
to express the short range correlation. 

As for the D-wave component, we have one Y2M func- 
tion in the spatial wave function. 

V'2Af(B£), {xi}) = exp(-ixBDx)|uxpy2M(wx) (9) 

Here, the global vector ux is defined as the linear com- 
bination of the vector x as {tx = X]f^^ ^UjX,;. The 
gaussian wave function with the function xB/jx = 
^ii^ ^ ^D^i ■ Xj" is able to optimize the role of the ten- 
sor interaction with small range gaussian functions in the 
d-wave part as the case of the TOSM. Since the spatial 
wave function now has a finite angular momentum, we 
have to make coupling of the spatial wave function with 
the spin wave function to get a desired total spin s. 

Nd 

\^D)^Y.^DMM^DA^^})xAW^})]sXti{n}))m 

We should know the role of the antisymmetrization 
operator A, which is defined as .4 = J^iLiPi^^i- 
sists of the particle exchange operator with the phase 
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Pi due to even or odd permutation, where the operator 
Vi makes the transformation of the relative coordinate 
xp = 7iX, the exponent of the correlated gaussian basis 
function and the generalized coordinate in the spherical 
harmonics takes the form after the antisymmetrization 
of the particles as Bp = 7T3T and up = Tu. The spin 
and isospin wave functions are subject to the antisym- 
metrization also. Thus, the anti-symmetrized wave func- 
tion contains quite a general expression of wave functions 
and the manipulation of various matrix elements is highly 
involved. We introduce the representation of the global 
vector for the basis wave functions in the same way as the 
Niigata group [7]]. In the stochastic variational method 
(SVM), the gaussian ranges B5 and Bd are generated 
randomly and choose the most suitable wave functions 
successively 

We obtain the wave function and energy in a few-body 
framework by taking minimization of the total energy 
"^■^T^^^y^ ~ ^^'^ ranges of the spatial wave functions 
Bs in Eq. ([8]) and Bpi in Eq. ([9]) and their amplitudes 
Cg and are chosen variationally. We may take also 
the intermediate spin and isospin quantum numbers in 
Eqs. (O, d?]) and PH]) with constraint that the final spin 
and isospin are s and t. Although there are several few- 
body methods, we work out numerical calculations using 
the SVM of the Niigata group 0]- We refer all the details 
of calculations in Ref. [3] . 

We show numerical results in the TOFM and compare 
with the rigorous few body calculation for A=2, 3 and 
4 systems and the TOSM calculation of Myo et al. for 
■^He The results for the deuteron shown in Table I 
and II indicate the importance of the tensor interaction. 
We work out numerical calculations with the TOFM us- 
ing the AV8' potential for A=3 and 4. We include only 
one relative coordinate with / = 2 angular momentum 
state ux = Xi, while the antisymmetrization A takes 
care of all the necessary permutations. We base our nu- 
merical calculation on the computer code of Varga and 
Suzuki 8] and add the tensor component in the varia- 
tional wave function. We obtain variationally suitable 
gaussian ranges B5 and Bp) and their amplitudes Cg 
in Eq. ([5]) and in Eq. (jlOp by diagonalization. In 
FiglU we show the energy convergence for as func- 
tions of the number TV of the basis wave functions. The 
energy components and the total energy are shown up 
to iV = Ns + Nd =150. We see a good convergence 
is achieved already around iV=50. All the energies are 
compared with the SVM calculation shown in the right 
hand side of the figure. We show in Figl^] the case of 
^He, which is obtained by using the 1=2 component with 
{tx = xi only. In this case the convergence is slow and 
we ought to go up to about 7V=300. The results of the 
rigorous few-body calculations for ^He are shown in the 
right hand side of this figure. Although the convergence 
is achieved for the central and the spin-orbit energies, the 
tensor interaction energy and also the kinetic energy are 
not yet converged. 

We show in Table I the total energy and various energy 




FIG. 1: Various energy components in the triton in the TOFM 
with AV8' as functions of the number A'^ = Ns + Nd of the 
basis wave functions. Those energy components are the ki- 
netic energy denoted by cross, the spin-orbit energy denoted 
by plus, the central interaction energy denoted by plus-cross, 
the tensor interaction energy denoted by open square and the 
total energy denoted by closed square. The corresponding en- 
ergy components obtained by rigorous calculations are shown 
in the right hand side. 
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FIG. 2: Various energy components in *He in the TOFM with 
AV8' as functions of the number of basis wave functions. The 
notation is the same as Fig[T] 



components as the kinetic energy, the central interaction 
energy, the tensor interaction energy and the spin-orbit 
interaction energy for the deuteron, -^11 and ^He. We 
do not include the Coulomb energy in this comparison. 
Although very small, we add other I = 2 components 
({tx = X2 and X3) in the calculation for this Table. The 
present result is compared with the rigorous SVM cal- 
culation by Suzuki et al. [7|. The results of various few 
body methods are essentially the same as the SVM re- 
sult 1]. As for *He, we also compare the results with the 
TOSM. 

We would like to discuss the comparison of the TOFM 
with the SVM for and ''He. The TOFM resuhs are 
compared almost perfectly with the SVM results for ■^H. 
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TABLE I: Various energy components with AV8'. 



Nucleus 


Energy 


Kinetic 


Central 


Tensor 


LS 


deuteron 


-2.23 


19.95 


-4.49 


-16.64 


-1.03 


^H(TOFM) 


-7.54 


46.67 


-21.98 


-30.47 


-1.95 


SVMf7] 


-7.76 


47.57 


-22.49 


-30.84 


-2.00 


*He(TOFM) 


-24.08 


95.53 


-54.61 


-60.95 


-4.05 


T0SM[4] 


-22.30 


90.50 


-55.71 


-54.55 


-2.53 


SVM[1J 


-25.92 


102.35 


-55.23 


-68.32 


-4.71 


TABLE IL Various 


energy components with G3RS. 


Nucleus 


Energy 


Kinetic 


Central 


Tensor 


LS 


deuteron 


-2.28 


16.48 


-7.29 


-11.46 


-0.00 


^H(TOFM) 


-7.61 


39.82 


-26.70 


-20.69 


-0.04 


SVM[7] 


-7.73 


40.24 


-26.80 


-21.13 


-0.03 


'*He(TOFM) 


-25.22 


84.83 


-66.21 


-43.66 


-0.17 


SVM[7| 


-26.05 


86.93 


-66.24 


-46.62 


-0.13 



For ■*He, the binding energy is -24.08 MeV, which is close 
to the total energy of -25.92MeV of the SVM. The cen- 
tral interaction energy and the spin-orbit energy are quite 
close to those of the SVM. On the other hand, the ten- 
sor interaction energy is about 8MeV smaller and the 
kinetic energy is about 7MeV smaller than the SVM re- 
sult. Hence, the tensor component and the kinetic energy 
are slightly underestimated. 

As for the comparison with the TOSM calculation for 
^He, the present calculation of the energy value is better 
than the TOSM result Large differences are found in 
the matrix elements of the kinetic energy and the tensor 
interaction. This difference should mean that the TOSM 
calculation can be improved by taking a more general 
UCOM correlation function. Some defect in the short 
range correlation may be seen also in the LS component, 
which is underestimated in the TOSM. 

In order to see the interaction dependence, we show the 
case of the G3RS bare NN-interaction in Table II [lo|. In 
this case, the tensor interaction is weaker than the case 



of the AV8' potential. Hence, now the agreement of the 
TOFM with the full calculation is impressive. We have 
calculated also the second 0+ state in ^He. The excitation 
energy comes out to be Ex = 18.5MeV, which is to be 
compared with the experimental data of 20.21MeV. The 
amount of the tensor energy for the is close to the 
one of '^H, which suggests the 3N-I-N structure. Excited 
states are described nicely in the TOFM, while there is 
a Faddeev model approach pTj . 

We have formulated a tensor-optimized few-body 
model (TOFM) in the spirit of the TOSM. The TOFM 
wave function contains the S-wave and D-wave compo- 
nents. The D-wave component contains only one Y2 com- 
ponent. Hence, the TOFM approximation makes the 
variational space much smaller that the one of the rigor- 
ous calculation. We have calculated A = 2,3 and 4 body 
systems in the TOFM and compared with the rigorous 
calculations. As for ^ = 3, the inclusion of single Y2 
component in the xi coordinate provides essentially the 
same results as the full model space calculation. This 
result indicates that one Y2 component is enough to take 
care of the tensor interaction in the three body system. 
As for ^ = 4, we again obtain good reproduction of the 
rigorous results, but the TOFM slightly underestimates 
the tensor interaction with single Y2 component. The dif- 
ference between the TOFM and the rigorous calculation 
comes essentially from the additional D-wave component 
in the other p-n pair. The inclusion of two Y2 functions 
in the variational space brings the total e nerg y very close 
to the rigorous result with AV8' for "'He [l^ . 

The present results indicate that nuclei like to have 
deuteron configurations and it is satisfactory to take the 
deuteron-like tensor correlations in the wave function. 
The present study is very encouraging to extend our 
study for nuclei with ^ > 5 in the TOFM and to de- 
scribe medium and heavy nuclei using the TOSM ap- 
proximation with better description of the short range 
correlation. 
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